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ABSTRACT 

Photometry data were collected from the literature and analyzed for supernovae 
that are thought to have a gamma-ray burst association. There are several gamma-ray 
bursts afterglow light curves that appear to have a supernova component. For these 
light curves, the supernova component was extracted and analyzed. A supernova light 
curve model was used to help determine the peak absolute magnitudes as well as es- 
timates for the kinetic energy, ejected mass and nickel mass in the explosion. The 
peak absolute magnitudes are, on average, brighter than those of similar supernovae 
(stripped-envelope supernovae) that do not have a gamma-ray burst association, but 
this can easily be due to a selection effect. However, the kinetic energies and ejected 
masses were found to be considerably higher, on average, than those of similar super- 
novae without a gamma-ray burst association. 

Subject headings: supernovae: general - gamma ray: bursts 



1. Introduction 



A significant effort was made, over a number of years, to find the origin of gamma-ray bursts 
(GRBs). The discovery of optical after glows in long-duration GR Bs allowed for the determination 
of the GRB's redshift (for a review, see Ivan Paradijs. et al.ll200Cl . and reference s therein). This led 



to the realization that lon g-duration GRBs are at cosmological distances (e.g.. Ivan Paradijs. et al 



19971 : iMetzger et al.lll997l ). In this paper, long-duration, cosmological GRBs will be referred to 
simply as GRBs. 

The light curve of the optical afterglow of GRB 980326 showed a rebrightening at around 
two weeks after the burst (jBloom et al.lll999l ). This rebrightening, or bump, in the light curve 
was attributed to an underlying supernova (SN). This was the first observational indication of an 
association between GRBs and SNe. 
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During April of 1998, an unusual gamma-ray burst, GRB 980425, was discovered w i th its 
gamma-ray energy several orders of magnitude lower than typical GRBs ([Soffitta et al 



1998 



Wooslev et al.lll999l ). Optical observations of the regi on led to the discove ry of an unusual super- 
nova, SN 1998bw, within the error radius of the GRB (jOalama et al.lll998l ). Spectra of SN 1998bw 
showed it to be a Type Ic SN with exceptional ly broad lines (Ic-BL). This SN was also unusually 
bright for a Type Ic SN (jRichardson et alJbood ). It was though t that GRB 980425 might simply be 



a normal GRB, but one wh ere the jet was viewed off axis (e.g. iNakamura 



1998 



loka Nakamura 



2001 



Yamazaki et al.l 120031 ). However, according to late-time radio observations (ISoderberg et al 



20041 ). even if it was viewed off-axis it would still be an unusual GRB. There may be more dim, 
peculiar GRBs, like GRB 980425, than observations would indicate; due to the fact that those at 
large redshifts are difficult to detect. 



The afterglow spectra of GRB 030329 were remarkably similar to that of SN 1998bw (jStanek et al 



20031 ). The SN associated with this GRB was then named SN 2003dh. Prom this, it was clear that 
at least some long-duration gamma-ray bursts are connected to supernovae. Since then, two other 
SNe Ic-BL have been spectroscopically confirmed to ha ve a GRB connection: SN 20031w/GRB 0312 03 



Modjaz et al 



2006 



(IMalesani et al.ll2004ll and SN 20 06ai/GRB 060218 (jSollerman et al 



Mirabal et al 



2006 



Mazzali et al. 2006b 



20061). Not all SNe Ic-BL have a convincing GRB connection 



as is the case for SN 1997ef (jlwamoto et al.ll2000l ) and SN 2002ap (jMazzali et al.ll2002l ). However, 
this does not rule out the possibility of an association with a weak (possibly off-axis) GRB. It also 
seems that there are some long-dur ation gamma,-ray b ursts with no apparent associated supernova; 
such as GRBs 060505 and 060614 (|Fvnbo et ahlbood ). 



Similar to the bump found in the optical afterglow light curve of GRB 980326 (mentioned 
above), there have been several o t her GRB afterglow light curves with just such a bump (e.g.. 

This provides us 



Reichart 



19991 : 



Galama et al. 



2000 



Sahu et al. 



2000 



Castro-Tirado et al 



2001 



with a number of possible GRB/SNe to consider, even without spe ctroscopic confirrnation. These 



Zeh et al 



20041 : 



Bloom 2004 



light curves are ty pically analyzed using a composite model (e.g 
Stanek et al.ll2005l ). discussed below in Section 5. Most of the light curves conside red here are of 
this type. For alternative explanations for this rebrightening see lFynbo et al.l (120041) and reference s 
therein. For a discussion on the possible progenitors of long-duration GRBs, see lFryer et al.l (|2007l ). 
The observed data are discussed in Section 2. Section 3 covers the analysis of the light curves while 
the results are presented in Section 4. The conclusions are discussed in Section 5. 



2. Data 

The photometry data were collected from the literature. R-band data were used for most of the 
GRBs due to their relatively high redshifts. V-band data were used for GRB 980425 {z = 0.00841) 
and GRB 060218 (z = 0.0331). The photometry references are given in Table 1. 



Distances and extinctions were taken into account in order to convert the apparent magnitudes 
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to absolute magnitudes. Luminosity distances were calculated from the redshifts using the following 
cosmology: Hq = 60 km s~^Mpc~^ , = 0.3 and Q\ = 0.7. The distance modulus, fi, for each 
GRB was calculated and is given in Table 1 along with the references. 



Foreground Galactic extinction was taken into account according to ISchlegel et al.l (|l998l ). 
These values were taken from NEE0 and converted from to Ay. For host galaxy extinction, 
the best estimates from the literature were used. In the case where the host galaxy extinction was 
estimated to be small or negligible (yet still unknown), a value of 0.05 mag was assigned. The 
extinctions are listed in Table 1, along with their references. 



Spectra of SI N" 1998bw were used to calculate K-corrections. 
Patat et al.1 (j200lh and obtained from the SUSPECT databas(i|. 



The spectra used were from 



3. Analysis 

Three of the light curves used in this study had no significant GRB afterglow component at the 
time of the SN: GRB 980425/SN 1998bw, GRB 031203/SN 20031w and GRB 060218/SN 2006aj. 
One of the light curves (GRB 030329/SN 2003dh) had no noticeable SN bump; that is to say that 
the SN component was maske d by a slowly decl ining afterglow. However, spectroscopy revealed 



a significant SN contribution (IStanek et al.ll2003l ). Most of the light curves had both a SN bump 



and a significant GRB afterglow component at the time of the SN. In order to analyze the SN by 
itself, the GRB afterglow component (and sometimes the host galaxy contribution) needed to be 
removed. If the light curve still included light from the host galaxy, and yet there was not sufficient 
late-time data to determine the host galaxy contribution, then that GRB was not included in the 
study. This was the case for GRB 991208. 

Once the host galaxy light had been accounted for, there were still two components. The GRB 
and the SN were treated independent of each other. While there is certainly some interaction be- 
tween the GRB and the SN, treating them as independent is a reasonable first order approximation. 
The early-time data, in the resulting light curve, were used to determine the contribution of the 
GRB afterglow. After day one, or day two, this was usually an unbroken power law. It appears as 
a straight line on a graph of R vs. log(t), as shown in Figure 1. The only exception is GRB 030329, 
where the GRB afterglow contribution was determined solely from the spectra. For all other GRBs, 
there needed to be sufficient early-time data to determine the GRB afterglow contribution. If this 
contribution could not sufficiently be determined from the data, then that GRB was not used in 
this study. This was the case for GRB 010921 and GRB 000911. 
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There were, however, a few exceptions (GRBs 020305, 020410 and 020903). The two dimmest 
SNe in the study are from the hght curves of GRB 020305 and GRB 020410. If the actual ORB 
contribution was larger than estimated for these SNe (shallower slope), then the SN contribution 
would be smaller making these SNe even dimmer than reported here. Even if the GRB contribution 
was negligible (steeper slope, including a downward break), then these two SNe would be brighter 
than reported here, but would still be the two dimmest SNe in the study. The light curve of 
GRB 020903 was not sufficient to get a good estimate of the GRB contribution. However, it was 
sufficient to reasonably determine that the GRB afterglow had a negligible contribution at the time 
of the SN's peak brightness. These three SNe are included because their peak absolute magnitudes 
are still relevant, even if the other information obtained from the light curve fitting remains highly 
uncertain (see Table 2). In general, the uncertainty in determining the GRB afterglow contribution 
is difficult to quantify and has not been included in the peak absolute-magnitude uncertainties 
given in Table 2. 



3.1. Model 



A SN light-curve model was used to help obtain accurate peak absolute magnitudes from the 
resulting light curves. It was also used in determining estimates of the kinetic energy, ejected mass 
and nickel mass for each SN. The mod el us ed here is a s emi-analytical model derived from two 
already existing models; lArnettI (|l982l ) and iJeffervl (|l999l ). At early times, the Arnett model is 
used; where the diffusion approximation is valid. At late times, the deposition of gamma-rays 
dominates the light curve, and the Jeffery model accounts for this. The basic assumptions are 
spherical symmetry; homologous expansion; radiation pressure dominance at early times; that ^^Ni 
exists and has a distribution that is somewhat peaked toward the center of the ejected matter. Also, 
optical opacity is assumed to be constant at early times and gamma- ray opacity is assumed t o be 
constant at late times. The combined model is described in detail by iRichardson et al.l (|2006l ). 



The model uses the SN's kinetic energy, ejected mass and nickel mass as parameters. After 
searching a grid of parameter values, a least-squares best fit was used to determine the most likely 
parameter values for each light curve. In order to improve the results, the ratio of kinetic energy to 
ejected mass was constrained. SN spectra were used to determine this ratio; however, spectra exist 
for only four of the GRB/SNe in the stu dy. These ratios are g iven in foe Mq^, wher e 1 foe = 10 ^^ 



erg. The values are 5.0 f or SN 1998bw jNakamura et al.l l200lh. 4.8 for SN 2003dh (iMazzah et al 



20031 ). 4.6 for SN 20031w (jMazzah et al.ll2006al ) and 1.0 for SN 2006aj (IMazzali et al 



2006J The 



average of these four values, 3.8, was used for the other GRB/SNe for which spectra were not 
available. Possible consequences of this approximation are discussed below (Section 4.2). 
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4. Results 

4.1. Absolute Magnitude Distribution 

The peak absolute magnitudes are listed in Table 2. The uncertainties shown were obtained by 
taking into account the uncertainties for each of the values in Table 1, as well as the observational 
uncertainties in the apparent magnitudes. Figure 2 shows a histogram of the absolute-magnitude 
distribution. This distribution has an average of My-pgafc = —19.2 it 0.2 and a sta ndard deviation of 



a = .7. Also shown in this figure is the distribution of stripped-envelope SNe, ([Richardson et al 



20061 . Fig. 2). Stripped-envelope SNe (SE SNe) are a combination of SNe lb, Ic and lib; and those 
shown here, do not have a GRB association. The main difference between these two distributions 
is that the GRB/SNe are, on average, brighter by 0.8 magnitudes. This is likely due to a selection 
effect. In order for most of the GRB/SNe to be detected, they have to be relatively bright compared 
to their GRB afterglow. Otherwise, it must be close enough to obtain a spectrum; as was the case 
with GRB 030329. Therefore, any relatively distant SN connected with a GRB that has a bright, 
or slowly declining, afterglow will not be detected; especially if the SN is relatively dim. This is 
why the dim GRB/SNe are likely to be undercounted. Notice that the two dim GRB/SNe fit well 
within the SE SN distribution. 

A graph of peak absolute magnitude versus distance modulus is shown in Figure 3. The two 
diagonal dashed lines represent the apparent magnitudes of 16 and 25. The horizontal dashed line 
represents the Type la ridge line, and is shown for comparison. The GRB/SNe can be compared 



to SE SNe (jRichardson et al.ll2006l ). included in this graph. Nearly all of the GRB/SNe were found 
to have an apparent magnitude dimmer than 16, but with a limiting magnitude of 25. This is in 
contrast to the SE SNe which were nearly all found to have an apparent magnitude brighter than 
16. Since this is related to distance, we see that it is rare to find any nearby GRB/SNe. Distant 
GRB/SNe are discovered because their associated GRBs are extremely bright. 

About half of the GRB/SNe in Figure 3 are near the Type la ridge line. The dimmest SN in 
the study, from GRB 020305, has a very large uncertainty. It is still worth including, however, due 
to the fact that it is firmly at the low end of the distribution even with the large uncertainty. 

When compared to other studies, the absolute magnitudes given here are, on av erage, somewhat 



bright er (after accounting for different cosmologies). For example, see T able 5 fromlSqderberg et al 



( 2006 ) and Figure 6 (which is similar to Figure 3 of this paper) from Ferrero et al. ( 20od ). This 



difference could possibly be due to the different methods for extracting the peak absolute magni- 
tudes. 



4.2. Light Curves 



The light curves are presented in Figure 4. Most of these light curves have a good range over 
time, considering the difficulty in separating the SN contribution from that of the GRB afterglow 
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(early times) and host galaxy (late times). 

The best estimates for kinetic energy, ejected mass and nickel mass are given in Table 2, along 
with rise times. Because spectra exists for only four of these GRB/SNe, E^/Mej values are only 
known for these four. The average value is used for the others. This is a reasonable first order 
estimate; however, a change in this value affects the individual and Mej values determined by 
the model. The trend is that increasing the E^/Mej ratio leads to an increase in both the E^ and 
Mej values obtained in the best model fits. This does not significantly affect the M^vi values and 
therefore does not affect My^peak- Als o notice, from Table 2, t hat My^peak and E^ are not correlated. 



This was the case for SE SNe as well (jRichardson et al.ll2006l ). Note that in the light-curve model, 



it is assumed that a substantial amount of ^^Ni is synthesized in the explosion. The decay of this 
isotope powers the peak of the light curve. Then, ^^Co (which is a product of the ^^Ni decay) itself 
decays and powers the tail of the light curve. 

The best fit model light curves are all shown on the same graph in Figure 5, with a common 
time of explosion. This shows, as expected, a general trend that the brighter SNe have broader 
light curves. This, however, is not always the case. GRB 030723 is very narrow compared to 
GRB 970228, yet they have similar peak magnitudes. 



5. Conclusions 

The average peak absolute magnitude was found to be higher for the GRB/SNe than for the 
SE SNe. However, in view of possible selection effects, the difference may not be significant. There 
were two GRB/SNe at the dim end of the distribution, but well within the distribution of SE SNe. 

Most of the SN data analyzed here were taken from GRB afterglow light curves. The GRB 
light was treated as being independent of the SN light and was removed so that the resulting light 
curve could be analyzed as a SN light curve. These resulting light curves fit quite well with a SN 
light curve model (Figure 4). 

GRB afterglow light curves that are suspected of having a SN component are usually analyzed 
by a different method. Usually the observed data are fit t o a composite model, where all of 



the components are represented; GRB, SN and host galaxy (jZeh et al.ll2004l . for example). The 
assumptions made are similar to those made with the separation method of this paper. The main 
difference between the two methods is in the way the SN is treated. The composite method starts 
with a light curve of SN 1998bw (adjusted for the GRBs redshift). Two of the five free parameters 
in the composite model are then used in the overall fit to describe the SN; one for the brightness 
and one for the width of the SN contribution. While the separation method uses SN 1998bw for 
K-corrections, a general SE SN model is then fit to the resulting light curve. The two methods 
are similar; however, the separation method used here allows for closer analysis of the the SN. 
Reasonable estimates of the kinetic energy, ejected mass, nickel mass an d rise times ar e found. 



However, the values of kinetic energy and ejected mass from other studies (jNomotd 120061 ) tend to 
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be larger by a factor of approximately two. 



The coincidence of the GRB date and the derived date of the SN explosion is another point 
of interest. In all but two cases the difference between the two dates was less than a week. For 
GRB 050525A, the SN explosion date is estimated to have occurred about 10 days before the GRB 
date, but by looking at the model fit and the observed data in Figure 4, we see that the model 
was not able to simultaneously reproduce the narrow peak and the bright, late-time data point. It 
appears that a more accurate explosion date for the SN would bring it closer to the GRB date. For 
GRB 011121, the SN explosion date is estimat ed to have occurre d about 9 days before the GRB 
date. Other studies have given similar results ([Bloom et al.ll2002l ). but the lack of pre-peak data 
for this SN makes the SN explosion date difficult to accurately pin down. Thus there is no clear 
evidence for real differences between the times of the GRBs and the SNe. 

The nickel masses determined here range from 0.05 to 0.99 Mm. These are similar to the values 
found for SE SNe with no GRB association ([Richardson et al.ll2006l ). However, the kinetic energy 
values determined here range from about one to 31 foe and the ejected mass values range from 
about one to six Mq. These values are, on average, considerably higher than those found for SE 
SNe with no GRB association. 
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Table 1. Data Used In Determining Absolute Magnitudes 



GRB/ 


Photometry 






Ref. 


Afi(Galactic)'' 


Ay {host) 


Ref. 


XRF 


Ref. 


(mag 


) 




(mag) 


(mag) 




970228 


1 


43.463 


± 


0.003 


2 


0.543 ± 0.087 


0.15 ± 0.15 


3 


980425 


4,5,6 


33.13 


± 


0.26 


7 


0.194 ± 0.031'= 


0.05 ± 0.05 


8 


990712 


9,10 


42.224 


± 


0.005 


11 


0.090 ± 0.014 


0.15 ± 0.1 


12 


011121 


13,14 


41.764 


± 


0.006 


14 


1.325 ± 0.212 


0.05 ± 0.05 


14 


020305 


15 


40.29 


± 


1.09 


15 


0.142 ± 0.023 


0.05 ± 0.05 


15 


020405 


16 


43.463 


± 


0.016 


16 


0.146 ± 0.023 


0.15 ± 0.15 


17 


020410 


18 


42.60 


± 


0.43 


18 


0.398 ± 0.064 


0.05 ± 0.05 


18 


020903 


19 


40.844 


± 


0.003 


19 


0.0935 ± 0.0150 


0.26 ± 0.26 


19 


030329 


20 


39.88 


± 


0.01 


21 


0.0675 ± 0.0108 


0.39 ± 0.15 


17 


030723 


22 


43.08 


± 


0.72 


23 


0.0886 ± 0.0142 


0.23 ± 0.23 


22 


031203 


24 


38.775 


± 


0.003 


25 


2.772 ± 0.444 


0.05 ± 0.05 


26 


041006 


27 


43.55 


± 


0.03 


27 


0.0607 ± 0.0097 


0.11 ± 0.11 


17 


050525A 


28 


43.10 


± 


0.07 


29 


0.2546 ± 0.0407 


0.12 ± 0.06 


30 


060218 


31 


36.15 


± 


0.05 


31 


0.471 ± 0.075'= 


0.13 ± 0.13 


31 



^Luminosity Distance {Hq 
ISchleeel et"al] l ll99d) 

'^Ai/ (Galactic) was used. 



: 60, Um = 0.3, Ha = 0.7) 



Refe renc es. — fl) |Galama et"al] l[2000h. (2) iBloom et all ll200ll') . ( 3) ICastander and Laiiibl 
lll999h. (4) iGalama et aTT 1998*'). (5) 'McKenzic fc Schac ferl (Il999l). (6) ISollerman et all ||2000|).^ 
(7) iKav et all lll998l) ,~(8) Nakamura ct al. (2001), (9) Sahu et al.l ||2000| '). (10) Fruchter et al, 
(200C), (11) NED. (12) IChristensen et al.l 112004). f 13 ) IBloom et al.l ll2002l). (14) iGarnavich et al' 
I 2OO3I). (15) iGorosabel et all ||2005(). fl6)lMasetti et all ||2003^. f 17 ) iKann et al.l ||2006|). (18) 



(251 



2003). (22) iFvnbo et al 



Levan et al] ||2005|). (19) iBersier et alt l l20od). (20) 



Prochaska et al 



Matheson et al.] ||2003|). 

Tominaga et al.l 120041) , (24) 



200. 



t004) . ( 23) ■ " . . - 1 

i ). (26) iMazzali et al.l I 2006e 



iDella Valle et al. I ll2006l ). (29) lFolev et al.l ||2005|) ,~(30) blustin et aLN2006l ). (31) ISollerman et all 



(27) 



21) iGreiner et atl 
Malesani et al.l l l2004ir 



Stanek et al. 



l 2005h. (28) 



(12006 
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Fig. 1. — The observed R-band light curve of GRB 990712 after the host galaxy light had been 
subtracted. The diagonal dashed line represents the contribution due to the GRB afterglow. In 
this particular case, the afterglow light can be described by the equation: Rag = 2.31log{t) + 21.22. 
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Fig. 2. — The peak absolute-magnitude distribution for GRB/SNe (solid bars) is shown with an 
average value of My^peak = —19.2 it 0.2 and a standard deviation of o" = 0.7. SE SNe are shown for 
comparison (striped bars). 
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Fig. 3. — Peak absolute-magnitude is plotted here versus distance modulus. The diagonal dashed 
lines are lines of constant apparent magnitude (16 and 25 mag.). The dashed horizontal line at 
My = —19.5, representing the SN la ridge line, is shown for comparison. The filled circles are 
GRB/SNe and the open squares are SE SNe. The associated GRB names are used to label a few 
key GRB/SNe. 
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Fig. 4. — All of the GRB/SN light curves in the study are plotted here with the best model fits. 
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Fig. 5. — All of the GRB/SN model light curves in the study are plotted here. The peak absolute 
magnitudes are given for each. 
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Table 2. Results 



GRB/ 


^V,peak 






Mjvi 




N 


XRF 


(mag) 


(foe) 


(Mq) 


(Mo) 


(days) 




970228 


-18.9 ±0.3 


23.2 


6.11 


0.54 


24 


4 


980425 


-19.4 ±0.3 


31.0 


6.22 


0.78 


23 


104 


990712 


-18.9 ±0.1 


5.32 


1.40 


0.20 


15 


5 


011121 


-19.6 ±0.3 


14.2 


3.73 


0.73 


21 


9 


020305 


-17.6 ±1.1 


3.61" 


0.95" 


0.05 


14" 


7 


020405 


-19.8 ±0.2 


11.1 


2.92 


0.72 


19 


5 


020410 


-18.0 ±0.4 


6.42" 


1.69" 


0.10 


16" 


3 


020903 


-19.5 ± 0.3 


11.9" 


3.13" 


0.60 


20" 


11 


030329 


-19.5 ± 0.2 


17.4 


3.63 


0.62 


20 


17 


030723 


-19.1 ± 0.8 


3.23 


0.85 


0.19 


13 


9 


031203 


-19.9 ±0.5 


21.0 


4.56 


0.99 


21 


10 


041006 


-19.9 ±0.2 


14.5 


3.81 


0.94 


21 


5 


050525A 


-18.9 ±0.2 


15.7 


4.14 


0.40 


21 


7 


060218 


-19.2 ±0.2 


0.89 


0.89 


0.30 


16 


42 



''These values are highly uncertain due to the difficulty in determining 
the contribution of the GRB afterglow. 



